Neutrino production in matter with time-dependent density or velocity 
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We show that neutrinos can be produced through standard electroweak interactions in matter 
with time-dependent density. 
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Electroweak interactions couple neutrinos to matter 
fermions, which have non-zero SU(2) charge. Both or- 
dinary matter and the nuclear matter of neutron stars 
carry net SU(2) charge. One can expect, therefore, that, 
if the matter density changes with time, neutrinos can be 
pair-produced, in analogy with the photon production by 
a time-dependent density of electric charge. We will show 
that this is, indeed, the case. 

Neutrino pair-creation in dense matter is similar to the 
pair-creation of fermions by a variable electromagnetic 
field H |, I or by an oscillating scalar field [|, |. The 
mechanism we describe is different from that discussed 
in 0, where the matter density gradients are essential. 
In contrast, we will assume the density is uniform and 
time-dependent. The discussion below is applicable to 
both ordinary matter and nuclear matter. For definite- 
ness, let us consider the latter. The weak interactions of 
neutrons n with neutrinos v are described by the follow- 
ing potential: 
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where Gp is the Fermi constant. In nuclear matter, the 
current j'^ = {GF/^/2){nj'^{l — 75)^), has a non-zero 
expectation value. For matter at rest, j^ = (p, 0,0, 0); 
here p = {GF/^/2)nn, and n„ is the number density of 
neutrons. For example, for matter densities pM ~ lO^*' — 
lO^^gcm"'^, inside a neutron star, the "reduced density" 
is p~ 10-3- 103 eV. 

If the matter density is time-dependent or if j changes 
direction, in the presence of interaction (|l|) , the neutrino 
field satisfies the Dirac equation with a time-dependent 
potential: 
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V- = 0. (2) 



Let us consider the case of a uniform time- 
dependent density p{t) moving with a four-velocity t;^ — 
{vo{t),0,0,V3{t)} along the third axis. Then j^ = pw^. 
We will assume that both p and w^ are constant at i < 0, 
and become time-dependent for i > 0. We now calcu- 
late the number of neutrinos produced by this variable 
nuclear current. 



Solutions of the Dirac equation can be written in terms 
of the cigcnmodes: 
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where the creation and annihilation operators satisfy the 
usual anticommutation relations. For the positive and 
negative energy modes we use the Ansatz: 
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and x± is an eigenspinor oi k ■ a with eigenvalue ±1. We 
will assume that vectors v and k are aligned because it 
simplifies the calculations; we do not expect the orthog- 
onal momenta to change the overall picture of particle 
production. 

Substituting this Ansatz into the Dirac equation, we 
get the equations for the mode functions fl ^{t), X = ±1. 
We introduce 

0(^i(i) = exp|-z/2jr*p±(C)dc}/SW, (7) 
which satisfies equation 
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where p± = p{vQ ± W3) for k aligned with the third axis. 
Equation (0) has two linearly independent solutions '*^ 



k,+ 



and 



kW 



for A = 1. Their complex conjugate counter- 



parts solve the same equation for A = — 1. In other words. 





FIG. 1: The occupation number of neutrino pairs, Afk'+, as FIG. 2: The number density of produced neutrino pairs, 



a function of time, for j° = 1 + 0.1sin(a;£) and j = 0, all in 
units of m. The upper curve corresponds to {u),k) = (10,4.3) 
(close to resonance), the middle curve to (1,-0.5), and the 
lower curve to (0.1, —0.5). 



In the presence of a matter background the charge con- 
jugation symmetry C is spontaneously broken. Equa- 
tion (H) relates the particles and antiparticles by a trans- 
formation that simultaneously changes the sign of p in 
equation (0). 

The initial state is specified at i < 0, when the den- 
sity is time-independent: /0±(t < 0) = 0. In a constant 
density background, the mode functions are plane waves, 
and we can choose the following initial conditions: 
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Here 
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are the eigenvalues of the Hamiltonian in matter with 
a constant density. In the limit p —> 0, the spinors U 
and V become the usual free field positive and negative 
frequency solutions. We normalize the spinors 



Ul ^{0)Uk,s{0) = Vl^JO)V-k.s{0) = Sr 



(12) 



Due to initial conditions, the orthogonality relation 

UlMV^kAO)^0 (13) 

is satisfied automatically. Since the time evolution is 
induced by a hermitian Hamiltonian, the normalization 
and orthogonality properties of the spinors are preserved. 
Condition ([l2) fixes the normalization constants: 
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To simplify notation we will drop the A label from the 
normalization constants, which are equal for A = ±1. 



J^k.+j £^s a function of fc at t = 10. The parameters chosen 
for this plot are j =1-1- 0.1 sin(a;t), j = 0, and lo — 1, all in 
units of m. 



The particle number rik^t) — {b\. g{t)bk,s{t)) is defined 
in terms of creation and annihilation operators. However, 
when bj, ^ and bk,s are time-dependent, the Hamiltonian 
does not, in general, remain diagonal at all times. This is 
a sign of particle production. To determine the number 
of particles, one has to diagonalize the Hamiltonian by a 
Bogoliubov transformation or, equivalently, project the 
states at time t from the basis of mode functions f{t) to 
that of plain waves. A detailed discussion and various 
applications of this standard technique can be found in 
Refs. Hill. 

We define the particle number in the basis of plain wave 
states in the background of constant density p± = p±(0). 
The positive and negative energy modes are given by 
eqs. (§,||J6|) with energies (|ll|) and normalization (|l4|). 
In this basis the time-dependent particle number is 



A/r.(i) = i^Li'' 
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where D^ ^ is the corresponding Bogoliubov coefficient: 



D^^ = 
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t/k,±^(0)y-k,±(t) 
= rnNl^ {^fL'\t) - /i:/^(t)4:i^(0) 
Therefore, the particle number is 
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In all formulas the neutrino momentum only enters in 
the combination kcg^±{t) — k±p±(t)/2. The momentum 
is shifted due to the effective potential generated by the 
matter background. 
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FIG. 3: The occupation number for neutrino pairs, J^k.+ ) 
as a function of fc at t = 175. The parameters are j^ — 
1 + 0.1 sin(ajt), j — and u = 5, all in units of m. The modes 
around k ~ 1.79 are in resonance. 



The number of particles at t = is equal to zero 
due to the initial conditions (HO) and (pT|). Further- 
more, equation (S) implies that at all times the number 
of particles A/jf^(i) is equal the number of antiparticles. 
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t(0)f/fe,±(i)|2, asitshouldbe. 



The total number density of neutrinos pairs is 



'{t) = 
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(27r)= 
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In the limit of an infinite volume the net amount of neu- 
trinos produced is n^{tcnd) — n^{to), which goes to zero 
if the matter density returns to its original state. How- 
ever, in reality the matter with time-dependent density 
occupies a finite volume, so that neutrinos can also diffuse 
away from the production region and break the coherence 
that was implicitly assumed in derivation of eq. (uTh . If 
the diffusion is fast, the density of neutrinos is smaller 
than the maximal equilibrium value (uTh and is nearly 
constant. In this regime, one can define a particle pro- 
duction rate niy{t) M. 

Particle production vanishes in the limit of a massless 
neutrino, see eq. (p7|). This is easy to understand in the 
current formalism. In the massless limit the right and 
left handed neutrino sectors decouple, and if the Hamil- 
tonian is diagonal at the initial time it remains so at all 
times. The Boguliubov transformation is trivial, and, 
hence, no particle production occurs. This is a conse- 
quence of the perturbation Sp{t) being spatially homoge- 
neous. The neutrinos are produced in pairs, and the two 
neutrinos have opposite momenta in the center of mass 
frame. Thus, the matrix elements of the form (P_gr^z/j), 
where V^ = 7'^, 1^1^ , vanish for massless neutrinos. Spa- 
tial perturbations 5p{x, t) lead to a particle production 
only if they are time-dependent in all Lorentz frames. For 
example, for a plane wave perturbation 5p ex cos(fc-x— wi) 
one can always make a Lorentz transformation such that 



6p' oc cos(fc'-x), and particle number remains constant in 
time. Here the phonons are virtual, and pair production 
is only possible if additional energy is pumped in the sys- 
tem, as is done in the production mechanism described 
in 1^. S Using the mode equations (||), we can com- 
pute the first and second derivatives of A/jf^ at i = 0+ . 
Let us consider time evolution of a two-component vec- 
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equation as 
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where 5 is a 2x2 matrix: 
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The expression ( |l7| ) can also be rewritten as 

Ml^{t)^m^Nl^T\t)Tnt), (21) 

where 
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Now we differentiate eq 
TT{0) = 0: 



k.± (0) 1 

ll) and use the identity 
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Nl^T^ (0) (S^T + TS)T{0) = 0; 



{S^'T + TS)S + {S^T + TS)]T{0) 
= 2m^Nl^ J^\Q)S^TSJ'{Q). 



Finally, 



AAi:±(0) = 0; A/-,^±(0) = -mXt±P±(0). 



(23) 



(24) 



Let us now consider a spatially uniform matter density 
that is changing periodically with time: 



jp(t) = p{t)Vf_„ Vf_, = const. 



(25) 



Figures 1-3 show the numerical solutions for the case 
j'^/m — 1 + O.lsin(wt), j = and various values of 
(u;,fc). The occupation number density ( p7| ) is a peri- 
odic function of time, as shown in Fig. 1. This means 
that neutrinos are created and absorbed through non- 
perturbative interactions with matter (c/. Ref. M). At 
all times Af^j_{t) remains smaller than 1, in accordance 
with the Pauli exclusion principle. The neutrino pair- 
production occurs in a spectral band centered around 
kcffito) — fc±p±(io)/2 — 0, as shown in Fig. 0. The value 
^cff (^0) — minimizes the real part of the frequency in 
eq. (0). In the adiabatic regime, where p ^ fc^jj -t- m^. 



the width of the resonance peak is F 



This clearly 



shows the non-perturbative nature of neutrino pair pro- 
duction in a coherent background. A perturbative cal- 
culation involving multiple phonons, each decaying into 
a neutrino pair, would predict a peak in the spectrum 
around fccff ~ w/2 with width F ~ GpUj^ for lu > m. 

For oscillations with p{t) = po + Spsm(ujt), Sp <ti po, 
equation (0) is a (complex) Mathieu equation with purely 
imaginary parameter q = iSp/uj and the ^-parameter 
that is also time dependent: A{k, t) — 4(m^ -t- {p± {t)/2 ± 
k)'^)/uj'^. Strictly speaking, this time dependence makes 
equation (S) different from the usual Mathieu equation, 
but the general behavior of its solutions is the same for 
Sp <^ p. Since the g-parameter is purely imaginary, all 
solutions are semi-periodic. There are several resonances, 
the first (and the widest) of which occurs when A = 
1, that is when the frequency and the momentum are 
related by iuj^ ^ m'^ + {k ± po/2)^. Fig. 3 shows this 
resonance for oj'-'^^ = 5m, at time t = 175m^^. The 
occupation number A/jf^ grows for those modes that are 
in resonance. 

There are additional resonant bands for higher values 
of A, that is for smaller oscillation frequencies. How- 
ever, since the ratio q/A oc u decreases with to, the 
higher resonances have diminishing width. Dense stars 
and other astrophysical objects are generally character- 
ized by small frequencies, so that neutrino production 
occurs off-resonance. 

The parameters in Fig. 1-3 were chosen so as to il- 
lustrate the general behavior of the solutions. We now 
examine the realistic values of these parameters in the 
neutron stars. 

Unless the lightest neutrino is much lighter than the 
0.1 eV scale inferred from the mass-squared difference 
measured by Super-Kamiokande, the basic oscillation fre- 
quencies in neutron stars have uj <$^ m, very far from the 
lowest - and widest - resonance. Of course, higher modes 
of oscillation may be excited, but they are usually as- 
sociated with significant dissipation. Let us consider the 
efficiency of neutrino production for p{t) = pQ+Spsm{ujt) 
with a small driving frequency lu <$^ m. Based on expan- 
sions of the Bogohubov coefficients in eqs. ( [l6| , p7|) , we 
estimate that the neutrino energy production per unit 
volume is e '--^ tti^ max[?TT,^, {Sp)'^]u). For a realistic case 
m ^ 6p '^ O.leV, one can estimate the neutrino produc- 
tion rate in a neutron star with radius i? w 10 km for the 



basic oscillation frequency uj '^ 2t:c/R: 
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This corresponds to emission of 10'^^ low-energy (sub-eV) 
neutrinos per second. For smaller amplitudes or larger 
frequencies, the production rates decrease or increase, 
respectively. 

This effect is, clearly, negligible when compared to the 
energy scale of a supernova explosion. However, it is con- 
ceivable that future observations of neutron stars may 
reveal subtle changes of the rotation speeds of cold oscil- 
lating neutron stars that can be compared with our pre- 
dictions. Neutrino production through the same mecha- 
nism can also occur in a neutron star binary system. In 
the latter case pair-production occurs because the cen- 
tripetal acceleration a ~ v j^ 0, and the current jp,{t) has 
a non-zero derivative. 

In summary, we have described pair creation of neutri- 
nos in a background with time-dependent matter density. 
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